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Abstract 

Using a collective- mode Monte Carlo method (the Wolff-Swendsen-Wang algorithm), 
we compute the spin-stiffness of the two-dimensional classical Heisenberg model. We show 
that it is the relevant physical quantity to investigate the behaviour of the model in the 
very low temperature range inaccessible to previous studies based on correlation length 
and susceptibility calculations. 

PACS numbers: 11.10Gh,11.10Lm,75.10Hk,75.30Fv 
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As well-known, the long-distance, low-energy physics of two-dimensional spin systems are 
expected to be obtained from a low-temperature perturbative expansion of a suitable Non 
Linear Sigma (NLc) model. However, the relevance of the low-temperature expansion beyond 
perturbation theory relies on the assumption of asymptotic freedom. Although this property 
is hardly questionable for symmetric 0(N)/0(N — 1) models with N > 2 (see however Jl|) no 
definite answer is known for more general models, in particular for Non-Symmetric models such 
as 0(3) <8> 0(2)/ 0(2) which are relevant for the study of frustrated Heisenberg spin systems 
0. Indeed, these models are suspected to be strongly affected by topological excitations [[| 
which is a possible source of failure of the low-temperature expansion. In order to have a non 
perturbative control of the low-temperature expansion, one can take advantage of Monte Carlo 
simulations. Up to now, calculations have been mainly concerned with correlation lengths 
and susceptibilities ||. Unfortunately, because of their exponential behaviour as a function of 
(3 = 1/kT and the computationaly accessible lattice sizes, studying the very low temperature 
regime is very demanding, or even impossible. The aim of this paper is to show that the 
relevant physical quantity allowing to reach this regime for accessible sizes is the spin-stiffness 
p s , a measure of the free energy increment under twisting of the boundary conditions ||, ||. 
In the following we shall restrict ourselves to the two-dimensional classical Heisenberg model; 
applications to more involved models will be presented in a forthcoming work. 

The Hamiltonian of the Heisenberg model is: 

H = -Jj2S l .S j (1) 

<ij> 

where < ij > denotes the summation over nearest-neighbours of a finite square lattice of size 
L. In (]]]), Si are three-component unit-length classical vectors and J is positive. Each site i of 
the lattice is indexed by two coordinates x» and 

We impose a twist in the x direction, by coupling the system with two walls of spins: 
S(x = 0) = Si, S(x = L) = S2, S2 being deduced from Si by a rotation of angle 6 around a 
direction e. The spin-stiffness p s is defined as: 



d*F(6) 

Ps(L) - 



(2) 



oe 2 

where F is the free energy. 

In terms of the spins it writes: 

p s (L) = L < £ (S^-S.e S r e) (x t -x 3 ) 2 > < ( £ (S.AS^.e ( Xl -x 3 ) f > (3) 

L <ij> 1 L <ij> 

where T is the temperature and Boltzmann averages are performed with two walls of parallel 
spins fixed at boundaries in the x direction. 
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The finite size behaviour of p s {L), when L is much larger than the lattice spacing a but 
much smaller than the correlation length £, has been calculated at one and two-loop order with 
use of the 0(3)/ 0(2) NLcr model [5,7]: 

1 £ 1 £ 
p s ln- + — In In- (4) 

Hs 2vr L 2it L K J 



where the common coefficient 1/2tt in front of the leading and sub- leading logarithmic terms is 
a universal number which is not modified by higher orders in the low-temperature expansion. 

The crucial point in measuring p s is that its predicted size dependence given by @ is all 
the more valid since L <C £. Therefore, in the very low temperature regime we can hope to 
test formula (|4]) by using a large range of relatively small lattice sizes. In contrast, measuring 
the temperature dependence of £ requires £ < L and therefore relatively high temperatures 
for accessible sizes[|J, a regime where the validity of the perturbation theory becomes less 
controlled. A most important point to notice is that at the very low temperatures considered 
here the physics of the model is entirely controlled by collective excitations - spin waves- and 
therefore we must take great care of these large-scale moves in any simulation of the model 
("beating" the critical slowing down). 

The purpose of this paper is to present a Monte Carlo study of the spin- stiffness for the 
finite two-dimensional classical Heisenberg model free of critical slowing down and then to 
investigate numerically prediction (£|). To summarize what have been obtained, our Monte 
Carlo calculations confirm the existence of a leading logarithmic contribution with the univer- 
sal amplitude 1/27T. In addition, an extra-contribution to the spin-stiffness consistent with the 
subleading term of has also been clearly identified. The Monte Carlo results presented have 
been obtained using the Wolff- Swendsen- Wang method || of updating large clusters of spins 
simultaneously. At the low temperatures considered here, using a collective Monte Carlo algo- 
rithm appeared to be essential to get well-converged values of the spin-stiffness. In particular, 
our preliminary attempts making use of a Monte Carlo algorithm based on local spin updates 
failed due to the severe critical slowing down. 

To our knowledge, we present the first unambiguous numerical calculation validating the pre- 
cise finite-size behavior of the spin-stiffness of the two-dimensional classical Heisenberg model. 
It should be noted that a similar calculation has been reported recently by Mon |J. However, 
we disagree both with the theoretical expression of the spin-stiffness used by the author and 
with the relevance of the local Monte Carlo scheme employed in his work. 

Results. The Wolff-Swendsen-Wang (WSW) algorithm has been implemented to simulate 
the Heisenberg model onaixL square lattice. In the y-direction periodic boundary conditions 
have been chosen. In the x-direction, fixed boundary conditions are to be used. Whereas in a 
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local Monte Carlo algorithm imposing fixed walls of spins pointing in some given direction is 
elementary, the situation is different when clusters of spins are built with WSW. To escape from 
this difficulty we have also chosen periodic boundary conditions in the x-direction. This intro- 
duces an error in the spin- stiffness exponentially small in InL. As expected, this contribution 
has been found to play no role in the following finite-size analysis. We have found that relatively 
moderate sizes L are in fact sufficient to validate formula (|4]). Lattices of sizes L=4,8,12,...32 
have been simulated. We have performed our simulations at four different temperatures: T/J= 
0.1, 0.15, 0.3, and 0.395. In each case we are at sufficiently low temperature to be in the regime 
of validity of formula (|) ( L « £). 

Figure 1 presents the complete set of results obtained for the spin-stiffness at different sizes 
and temperatures. At the scale of Figure 1, all curves appear to be very rapidly linear as 
a function of InL. In order to determine accurately the corresponding slope a closer look is 
necessary. Figure 2 presents a blow up of data of Figure 1 for the lowest (upper figure) and 
highest (lower figure) temperatures treated, T/J = 0.1 and T/J = 0.395, respectively. A first 
point to notice is that a very high accuracy on our data has been achieved. Such a level of 
accuracy is absolutely necessary to put into evidence the linear regime of the spin-stiffness as 
well as to get a truly converged estimate of the slope. We emphasize that only when resorting 
to a collective Monte Carlo scheme we have been able to fulfill both requirements. In our first 
attempts to use a local Monte Carlo scheme we observed a systematic and uncontrolled long- 
term drift of the estimates of the statistical mean values. A first important remark concerning 
Figure 2 is how fast we enter the linear regime: at all temperatures considered it is reached 
at L ~ 16. By using data for L = 16, 20, 24, 28, and 32 an estimate of the slope can be 
extracted, we get: -0.162(4),-0.166(5),-0.171(5), and -0.184(7) at T/J=0.1, 0.15, 0.3, and 0.395, 
respectively. At the very low temperature T/J = 0.1 we recover within statistical fluctuations 
the theoretical result 1/27T = 0.1592... predicted by formula (f|). [|l(| At higher temperatures 
non-negligible higher-order contributions in the spin-stiffness show up. To put this on a more 
quantitative basis, we have performed a fit of the data using the full expression ([|). The 
resulting curve is represented by a solid line in Figure 2. The only free parameter entering the 
fit is the correlation length £, the arbitrary reference value for the spin-stiffness being chosen 
so as to reproduce exactly the last data (L = 32). The dashed line is the linear curve obtained 
when resorting to the leading logarithmic behaviour (no In In corrections, no renormalization of 
the 1/27T slope) using the very same correlation length as determined in the fit. At T/J = 0.1, 
both curves almost coincide in the linear regime, illustrating the correctness of the leading 
log prediction and the smallness of the higher-order corrections at this temperature. At the 
higher temperatures considered, we clearly see the necessity of going beyond leading order. In 
addition, it is striking to see how good representation (|4]) is in reproducing our Monte Carlo 
data. Of course, at the accuracy determined by statistical fluctuations it is not realistic to 
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hope to resolve the precise analytical In In behavior of the second-order theoretical expression. 
However, our data are perfectly consistent with the "renormalized slope" predicted by (^), 
s * = dp s /d\nL = -1/2tt(1 + l/ln(f/L)). 

In Figure 3 we have plotted the correlation length £ issued from the fit using formula We 
also present the curve obtained from the formula proposed by S.H. Shenker and J. Tobochnik 
ITf (obtained by matching high- and low-temperature calculations): 



exp27rJ/T . . 

It is very satisfactory to see that our rough estimates of £ are in good agreement with this 
completely independent calculation of the correlation length. 
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Figure Captions 



Fig.l Spin-stiffness for different sizes and temperatures. Statistical fluctuations smaller 
than the size of crosses. 

Fig.2 Blow up of Fig.l for T/J = 0.1 and T/J = 0.395. The solid line is the best fit using 
Eq.|], the dashed line the first-order prediction (no renormalization of the slope). 

Fig.3 Correlation length £. The solid line is obtained from Eq.|, the values indicated by 
crosses from the fit of our data using Eq.|[ 
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